February 7, 2008 



Chaos and Non- Archimedean metric in the 
Bernoulli Map 

Jesus San-Martin 1 * Oscar Sotolongo-Costa 2 ' 3 
1- Dpto. Matematica Aplicada, Esc. Universitaria de Estadistica, 
Universidad Complutense. 
Av. Pta. de Hierro S/N. 28040 Madrid. 
2- LCTDI, UNED, a.p. 60141, 28080 Madrid, Spain. 
3- Fac. de Fisica, Universidad de La Habana, Habana 10400, Cuba. 



Abstract 

Ultrametric concepts are applied to the Bernoulli Map, showing the ade- 
quatencss of the non- Archimedean metrics to describe in a simple and direct 
way the chaotic properties of this map. Lyapunov exponent and Kolmogorov 
entropy appear to find a simpler explanation. A p-adic time emerges as a natural 
consequence of the ultrametric properties of the map. 
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I. Introduction 



After the work by Mezard et al.,0 ultrametricity has triggered the interest in 
a wide range of physical phenomena, due to its applications in different topics: 
spin glasses, mean field theory, turbulence, nuclear physics . Also optimization 
theory, evolution, taxonomy, protein folding benefits from it ( for an excel- 
lent review seeau) . Wherever a hierarchical concept appears, non- Archimedean 
analysis is an adequate tool to study the problem. 

Ultrametricity is a promising tool in the theory of branching processes, 
which alsaiias revealed its possibilities in the study of self-organized critical 
processes. H'H'El It seems possible to find simpler tools to describe the geometry 
of these processes. Here, we expose the advantages of a hierarchical represen- 
tation in the case of the Bernoulli shift. This will permit, using simple geo- 
metric considerations, to determine the magnitudes governing the system, and 
the advantages of a p-adic metric will be stressed over the euclidean one. The 
ultrametric distance will be shown to be consistent with the characteristic be- 
havior of this chaotic unidimensional map. In ultrametric spaces, concepts such 
as exponential separation of neighboring trajectories, and characteristic param- 
eters (Lyapunov exponents and Kolmogorov entropy) seem to find a simpler 
explanation than with the euclidean metric. 

As an example, where euclidean metric is not very adequate, let us consider 
the Baker's map.Q The interval [0, l]x[0, 1] is mapped to [0, l]x[0, 1]. There- 
fore, the distance between two points can't be larger than the distance between 
two opposite corners in [0, l]x[0, 1]. Nonetheless, the Baker's map has got a 
Lyapunov exponent bigger than one. Then, the distance between neighboring 
points grows exponentially in a finite region of the phase space. In the euclidean 
space we would have to define the distance in this case as the euclidean length 
of the shortest path lying entirely within the region that has suffered the defor- 
mation.tl As any nontrivialj-norm is equivalent to the euclidean or any of the 
p-adics (Ostrwski's theoremcl), it would be convenient to measure the distance 
between points in the Baker's map with a p-adic metric. 

An ultrametric space is a space endowed with an ultrametric distance, de- 
fined as a distance satisfying the inequality 

d(A, C) < Max{d(A, B) + d(B, C)} (1) 

( A, B and C are points of this ultrametric space), instead of the usual 
triangular inequality, characteristic of euclidean geometry 

d(A, C) < d(A, B) + d(B, C) (2) 

A metric space E is a space for which a distance function d(x, y) is defined 
for any pair of elements (x, y) belonging to E. 
A norm satisfying 
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\\x + y\\<max{\\x\\,\\y\\}, (3) 

is called a non- Archimedean metric, because equation 3 implies that 

\\x + x\\ < \\x\\ (4) 

holds, and equation 4 does not satisfy the Archimedes principle: 

\\x + x\\>\\x\\. (5) 

A metric is called non- Archimedean or ultramctric, if 

d(x, z) < max {d(x , y) , d(y , z)} . (6) 
A non- Archimedean norm induces a non- Archimedean metric: 

d(x, z) = \\x — z\\ < max{d(x, y), d(y, z)}. (7) 

Equation 7 implies a lot of surprising facts, e.g., that all triangles are isosceles 
or equilateral and every point inside a ball is itself at the center of the ball, 
furthermore the diameter of the ball is equal to its radius. 

An example of ultrametric distance is given by p-adic distance, defined as 

d p (x,y) = \\x-y\\ p (8) 
where the notation defines the p-adic absolute value: 

M\ p =p-\ (9) 

where p is a fixed prime number , x ^ is any integer, and r is the highest 
power of p dividing x.Two numbers are p-adically closer as long as r is higher, 
such that p r divides \\x — y\\. Amazingly, for p = 5, the result is that 135 is 
closer to 10 than 35. 

Any positive or negative integer can be represented by a sum 



E^ 1 ' ( 10 ) 

where 



x 

i=0 



0<Oi<p-l. (11) 

If negative exponents are considered in the sum, rational numbers can also 
be represented. Such a representation is unique. The set of all sums Q p is the 
field of p-adic numbers, and contains the field of rational numbers Q but is 
different from it. 
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II. Lyapunov exponent and Kolmogorov entropy 

With the above description the p-adic numbers have a hierarchical structure, 
whose natural representation, is a tree. Let us now use this description to work 
with the Bernoulli map (Seed). 

In this way, the numbers can be represented as a set of points in a straight line 
or by a hierarchical structure, depending on the definition of distance (euclidean 
or Archimedean). 

Let us represent the initial value (state) to be mapped into the unit interval 

by the sequence 0,ai, a^, with a,i — or 1 to denote the initial value in 

binary notation. 

It is possible to reorder these sequences as a hierarchical tree. To get it, 
let us do the following process to represent the result of the application of the 
Bernoulli map: 

We begin at an arbitrary point. We read, consecutively, the values of 

, from i=l to N, of the sequence 0, a% a^. When a, takes the value we 

move to the left, and the same distance down. When asj takes the value 1 we do 
the same, but moving on the right. The result is 2 N branches of a hierarchical 

tree. Any branch represents univocally a possible sequence 0,ai ajv. Thus, 

for instance, the sequence 0,0110 represent: left, right, right, left. Later, we will 
make N — > oo to recover the Bernoulli map. 

The distance d(xi,Xj) between two branches (sequences) x 1 •, x j m this tree 
is given by 

d(x i ,x j ) = { Q,i = i ' (12) 

where m is the number of levels one must move up the tree to find a common 
branch linking Xi and Xj , and N is the number of levels (the length of the 
sequence). This is equivalent to 

d(x i ,x j ) = { l/lf 3 , (13) 

where h is the position of the last block in which at (i = l,....,h) are 
common to the two sequences Xi, Xj. It means that the numbers Xi and Xj are 
close up to the h th binary place. This distance is an ultrametric one. 

Now we consider the map, given by /(0, a\ ajy) = 0, a-i a^a , such that 

aj, i=l,...,N moves one position to the left, and a new figure a is born, given 
the accuracy of the descripition. In the limit when N— > oo the Bernoulli shift 
will be retrieved. 

To calculate the Lyapunov exponent it is necessary to know how neighbor- 
ing points xq + e and xq evolve during the Bernoulli map. Let e be equal to 

2- h (l + 2~ Sl + 2~ &2 + ) > 2~ N , then the first different position between 

xq = 0, ai..ah—i a h—Q'N an d xq + e is ah- Then, it is necessary to move up the 
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tree N - h +1 levels from the bottom line to find the common branch in the 
position a/j_i (obviously, the last common figure between xq and xq + e). Then 



d(xo + e lX o) = 2-' ,+1 (14) 

and 

d(/"(x + e),r(xo)) = 2- /i+1 +" (15) 

because the iteration /" moves away the common branch n positions from 
the bottom level. 

To calculate the Lyapunov exponent it is necessary to express the exponential 
growth of the distance between two neighboring points: 

lim lim2 A "e= lim lim d(f n (x Q + e) , f n (x )) (16) 

n — >oo e — >0 n — >oo e — >0 

Since the base for measuring the p-adic distance in our space is the number 
2, in the preceding equation we have expressed the exponential growth with 2 A ™ 
instead of e Xn . 

Replacing e and d(f n (xo + e), f n (xo)) in the preceding equation we obtain 



lim lim 2~' l (l + 2~ Sl +2^ 2 + )2 A " = lim lim 2~' l+1+n (17) 

n^oo h—>oo n—>oo h—>CG 

from 17 it can be easily observed that A = 1. 

Since the Lyapunov exponent in Bernoulli map is In 2,0 we recover this results 
with p-adic metric, since 2 1 = e ln2 . That means that each unit time interval 
implies a new doubling of branches in each node of the hierarchical tree. Then, 
once a unit time interval has elapsed, the number of levels one must move up 
the tree to find a common branch increases in one. This result will be crucial 
to understand how the information is lost in the course of time. 

The former explanation has stressed the simplicity of the new expression 2 1 , 
in comparison with e ln2 , and its interpretation. 

In unidimensional maps, as the one considered here, the Kolmogorov entropy 
coincides with the Lyapunov exponent .0 The Kolmogorov entropy expression 
is: 



k = - lim lim — V" p H i n \g 2 Ph-i n , (18) 

n-too T-tO TIT L — ' 



H. 



where pi t , n is the probability to reach the i n state of the system in the 
phase space following a given path i\i%...i n . It can be seen that in our case 
this probability only depends on the final state i n because for each state there 
is just one path, i.e., that given by the sequence i\ii...i n . Besides, the number 
of states in the n th level is 2™ , and r is the time elapsed to pass from one 
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state to a successive one. The probability to occupy one of the 2™ states is 
Pn = Pi 1 i 2 ...i n = and it results 

2 

k = lim lim . (19) 

But the distance between two successive states of the n th level is 2 1_ ™, 
because they are common until the (n — l) th level. Since the speed v to pass 
from one sequence to the next is constant in Bernoulli map, i.e., v = 2 = 1 
the time r elapsed between these two successives states is r = 2 1_n . As expected 
k = 1, coinciding with the Lyapunov exponent. 

The Kolmogorov entropy measures the loss of information in the process. 
From our representation this loss of information can be easily seen, since the 
process of separation of trajectories is such that for any step the increase of the 
distance between two points duplicates the number of branches through which 
this increment can be reached. We are loosing information because we don't 
know exactly the way we are separating two states. 

On the other hand, we can see that in the ultrametric space the natural time 
of the system is also ultrametric. The time of transition between two sequences 
Xi,Xj satisfies the same expression 12 as the distance between Xi,Xj. 

Besides, subsequent behavior of two states that separate in a given point in 
the ultrametric space depends from the point in which this separation occurs, 
revealing an aging effect. This effect will be treated in future works. 



III. Conclusions 



It was verified that the Bernoulli map leads to a hierarchical structure in 
the p-adic metric. With the ultrametric distance the Lyapunov exponent and 
the Kolmogorov entropy acquire a simpler expression and a direct geometric 
interpretation is supplied by the hierarchical structure. The p-adic metric seems 
to be the natural metric of this map. The hierarchical structure generates p-adic 
properties for the temporal evolution. 
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